Two-level correlation function of critical random-matrix ensembles 
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The two-level correlation function Rd,p(s) of d-dimensional disordered models (d — 1, 2, and 3) with 
long-range random-hopping amplitudes is investigated numerically at criticality. We focus on models with 
orthogonal (J3 — 1) or unitary (J3 — 2) symmetry in the strong (b d <C 1) coupling regime, where the parameter 
b~ d plays the role of the coupling constant of the model. It is found that Rd,/3(s) is of the form Rd,/}(s) = 
1+6(b)— F/}(s^/b where Fi(x) = eiic(ad,px) andf^x) = exp(— ad, px 2 ), with ad,p being a numerical 
coefficient depending on the dimensionality and the universality class. Finally, the level number variance and 
the spectral compressibility are also considerded. 

PACS numbers: 71.30.+h, 72.15.Rn, 71.55.jv. 05.40.-a 



I. INTRODUCTION 

Random-matrix theories are largely and successfully ap- 
plied in the theoretical description of complex nuclei 1 ' 2 , gauge 
field theories* 3 * 4 * 5 , mesoscopic systems, 6 - 7 - 8 - 9 and random sur- 
faces in the field of quantum gravity. 10 Their advantage is that 
it is possible to represent the Hamiltonian of the correspond- 
ing system by a large Hermitian matrix acting on a finite- 
dimensional Hilbert space (if one disregards the continuous 
part of the spectrum). 

The random matrices fall into one of three universality 
classes, named orthogonal (j3 = 1), unitary (J3 = 2), and 
symplectic (f3 = 4), depending on the global symmetry prop- 
erties of the Hamiltonian they represent. 11 The symmetry pa- 
rameter (3 is the number of independent real components that 
characterizes a matrix element of the Hamiltonian. A system 
belongs to the orthogonal class if it has both time-reversal and 
spin-rotation symmetries, to the unitary class if time-reversal 
symmetry is broken, and to the symplectic class if the system 
has time-reversal symmetry but spin-rotation is broken. The 
relevant terms in the Hamiltonian are a coupling to an applied 
magnetic field, which breaks time-reversal symmetry, and the 
spin-orbit interaction, which breaks spin-rotation symmetry. 

One of the most relevant applications of the random ma- 
trix ensembles is to the study of critical phenomena, particu- 
larly to the special case of critical statistics which is found at 
the Anderson metal-insulator transition (MIT) in disordered 
systems, 12 . 13 . 14 i 15 . 16 . 17 i 18 . 19 . 20 i 21 . 22 . 23 i 24 . 25 At the critical point, 
the states acquire the property of multifractality, which marks 
a qualitative difference from the extended states in a metal 
and localized states in an insulator. These critical states corre- 
spond to critical-level statistics. Although several ensembles 
of nonconventional random matrices have been suggested to 
describe this statistic ^^J^^^i^k we wish to emphasize 
the power-law random-banded matrix model (PRBM)- 3 - 3 . for 
which the multifractality of eigenstates has been rigorously 
proven. 33 - 34 This model is characterized by a variance of their 
off-diagonal matrix elements, which decay as a power law 
with increasing distance from the diagonal. It should be men- 
tioned that all these models are of a one-dimensional nature. 

Energy-level correlations provide general tools for the sta- 
tistical description of disordered systems, helping in our un- 
derstanding of the localization transition. An important statis- 



tical measure of spectral correlations is the two-level correla- 
tion function (TLCF) of the density of states (DOS), which 
measures the correlations of the DOS at two different en- 
ergies. This function has been derived analytically for the 
PRBM model in the two limiting cases of weak and strong dis- 
order by mapping the corresponding Hamiltonian onto an ef- 
fective a model of a one-dimensional (ID) nature 35 and using 
renormalization-group methodsp^i respectively. We stress 
that, unlike the ID PRBM model, it has not until now been 
possible to analytically solve the most interesting, from an ex- 
perimental point of view, disordered models with long-range 
transfer terms in d = 2 and 3. Thus, explicit results for the 
TLCF in two-dimensional (2D) and three-dimensional (3D) 
models are still lacking, and finding this function is essential 
in order to fully understand the MITs. For this reason, we 
addressed the problem using numerical calculations. An im- 
portant and closely related quantity, the level-number variance 
(LNV), will be also considered. 

In this work we numerically calculate the TLCF and the 
LNV of critical d-dimensional random matrix ensembles with 
long-range off-diagonal elements and orthogonal or unitary 
symmetry. Since MITs genetically take place at strong dis- 
order (conventional Anderson transition, quantum Hall tran- 
sition, transition in d = 2 for electrons with strong spin-orbit 
coupling, etc.), we will restrict ourselves to the study of both 
quantities in this regime. In the ID case our results for the 
TLCF are in good agreement, except for the numerical co- 
efficient, with existing analytical estimates. 35 In addition, we 
propose expressions for the TLCF in the d — 2 and 3 cases. 
Apart from the importance of these findings from a general 
point of view, they may be relevant for several real physical 
systems (see Sec. Hill. 

The paper begins by first describing the model and the 
methods used for the calculations in Sec. |II] The results for 
the TLCF and the LNV in models with orthogonal or unitary 
symmetry are presented in Sees. llllland HVl respectively. Fi- 
nally, Sec. IVlsumm arizes our findings. 



II. MODEL AND METHODS 

In order to fully represent the mesoscopic systems we intro- 
duce an explicit dependence on dimensionality d in the widely 
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studied PRBM ensemble, 33 ! 34 ! 35 . 38 ! 39 ! 40 ! 41 . 42 ! 43 . 44 . 45 ! 46 . 47 . 48 ! 49 . 50 
Thus, we consider a generalization to d dimensions of this 
ensemble. The corresponding Hamiltonian, which describes 
noninteracting electrons on a disordered <i-dimensional square 
lattice with random long-range hopping, is represented by ran- 
dom Hermitian L d x L d matrices Ti (real for (3 = 1 or complex 
for f3 — 2), whose entries are randomly drawn from a normal 
distribution with zero mean, (7Uj) — 0, and a variance that 
depends on the distance between the lattice sites r. 
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in which standard Gaussian ensemble normalization is used. 51 
Using field-theoretical methods, 26 . 33 . 34 . 36 ! 38 ! 39 ! 52 ! 53 the 
PRBM model was shown to undergo a sharp transition at 
a = d from localized states for a > d to delocalized states 
for a < d. This transition shows all the key features of the 
Anderson MIT, such as multifractality of the eigenfunctions 
and nontrivial spectral compressibility at criticality. In what 
follows, we focus on the critical value a = d. 

The parameter b d in Eq. Q is an effective bandwidth that 
serves as a continuous control parameter over a whole line 
of criticality, i.e, for an exponent equal to d in the hopping 
elements Hi 



b d w& Furthermore, it determines the critical 



dimensionless conductance in the same way as the dimension- 
ality labels the different Anderson transitions. Each regime is 
characterized by its respective coupling strength, which de- 
pends on the ratio ({\Hu \ 2 )/{\Ti l0 1 2 )) 1/2 oc b- d between 
diagonal disorder and the off-diagonal transition matrix ele- 
ments of the Hamiltoniani 5 ^ 

Many real systems of interest can be described by Hamilto- 
nians Q. Among such systems are optical phonons in disor- 
dered dielectric materials coupled by electric dipole forces, 55 
excitations in two-level systems in glasses interacting via 
elastic strain, 56 magnetic impurities in metals coupled by an 
r -3 Ruderman-Kittel-Kasuya-Yodida interaction, 57 and im- 
purity quasiparticle states in two-dimensional disordered d- 
wave superconductors. 58 It also describes a particle moving 
fast through a lattice of Coulomb scatterers with power-law 
singularity, 53 the dynamics of two interacting particles in a 
ID random potential, 59 , and a quantum chaotic billiard with a 
nonanalytic boundary^ 

The TLCF is defined in the usual way 

Rd,p(u) = 7A2 ^ + u/2)v(e-u/2)), (2) 

where v(e) = L- d TrS(e - H) is the fluctuating DOS and ( ) 
denotes averaging over disorder realizations. At the critical 
point Rd^oj) acquires a scale-invariant form, if considered as 
a function of s = u>/A, the frequency normalized to the mean 
level spacing A = 1 / L d (v{e)) tiSu&llln the case of constant 
average DOS, Rd,p(s — lo/ A) can be simply rewritten as 



0) 



where p(n; s) is the distribution of distances s„ between n 
other energy levels and the S(s) function describes the self- 
correlation of the levels^ 



The strong disorder limit (b < 1) of the ID PRBM model 
can be studied using the renormalization-group method of 
Refs. and |^ For orthogonal symmetry, the following 
result is obtained for the TLCF at the center of the spectral 
bandit 



(4) 



i?i,i(s) = 1 + 5(s) - erfc ( ai (1 — 



where erfc(x) = (2/^/tt) ^ exp(—t 2 )dt is the complemen- 
tary error function and a\ t \ = 1/t/tt, whereas for unitary 
symmetry 



RiA s ) = 1 + S ( s ) ~ ex P 



-01,2 



b 2 



(5) 



with ai,2 = 2/n. Note that for small s, Ri^(s) behave as s@ 
thus reflecting the levels repulsion effect, whereas they tend 
asymptotically to 1 at large values of s. 

Another important quantity, which characterizes fluctua- 
tions in the level density on larger scales than the mean spac- 
ing, is the variance 



(6) 



of the number of levels in an energy window that contains 
1 -C (n) <C N on average. This variance is given in terms of 



the TLCF by E^((n)) - !% } ({n) - \s\)R d , p (s) ds. The 
number variance is a statistical quantity that provides a 
quantitative measure of the long-range rigidity of the energy 
spectrum. For the Poisson distribution, the levels are uncor- 
rected and there are large level-number fluctuations, leading 
to a linear variance g((n)) = (n). On the other hand, 
the level correlations in the Wigner-Dyson statistics make the 
spectrum more rigid and the number variance grows only 
logarithmically T? d p((n)) ~ ln(n). However, in the criti- 
cal regime, the variance has been conjectured to be Poisson- 



J d.3 



((n)) - x(n) , 



(7) 



where the level compressibility x is another important pa- 
rameter for characterizing the MIT and which takes values 
< x ^ 1> with zero referring to delocalized states and unity 
the localized states. 

For the computation of Rd,p(s) and p((ri)), we unfold 
the spectrum in each case to a constant density and rescale 
it so as to have the mean spacing equal to unity. Then we 
calculate p(n; s) and use Eqs. (0 and 0, respectively. The 
system sizes range between L — 500 and 6000 in ID, L = 20 
and 100 in 2D, and between 8 and 14 in 3D, whereas b d ranges 
in the interval 0.02 < b d < 0.12. We consider a small energy 
window, containing about 10% of the states around the center 
of the band. The number of random realizations is such that 
the number of critical levels included for each L is roughly 
1.2 x 10 6 . In order to reduce edge effects, periodic boundary 
conditions are included. 
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FIG. 1: Ri,i(s) for the ID model Q with orthogonal symmetry 
(/3 = 1) as a function of the rescaled variable s/b at different 6 
values for several system sizes: b — 0.02, L = 4000 (circles), 
b = 0.05, L = 6000 (squares), 6 = 0.08, L = 1000 (diamonds), 
b = 0.1, L = 4000 (up triangles), and b = 0.12, L = 1000 (left 
triangles). The solid and dashed lines represent the renormalization- 
group estimate [Eq. J4j], with the fitting parameter ai ; i = 0.502 and 
the predicted ai,i = 1/^/tt, respectively. The inset shows the same 
data on the scale s. 
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FIG. 2: Rd,i(s) for the 2D (upper curve) and 3D (lower curve) 
models with orthogonal symmetry (3 — 1 as a function of the 
rescaled variable s/b d at different b d values for several system sizes: 
b 2 = 0.02, L 2 = 40 2 (circles), b 2 = 0.05, L 2 = 20 2 (squares), 
b 2 = 0.08, L 2 = 60 2 (diamonds), b 2 = 0.1, L 2 = 100 2 (up trian- 
gles), b 2 = 0.12, L 2 = 30 2 (left triangles), b 3 = 0.02, L 3 = 12 3 
(circles), 6 3 = 0.05, L 3 = 14 3 (squares), b 3 = 0.08, L 3 = 12 3 (dia- 
monds), b 3 = 0.1, L 3 = 8 3 (up triangles), and b 3 = 0.12, L 3 = 10 3 
(left triangles). The solid lines are fits to Eq. with the fitting 
parameters d2,i = 0.308 and as t i — 0.208. 



III. TWO-LEVEL CORRELATION FUNCTION 

In this section, we numerically compute the TLCF Rd,@{s) 
of Hamiltonians Q with the orthogonal or unitary symmetry 
for different values of the inverse coupling constant b d ^ 1 
and various system sizes. We also compare our results with 
the analytical estimates of Ref. 35 for the ID model. 



A. Orthogonal symmetry 

Let us first check the renormalization-group result, Eq. 
@, corresponding to the ID Hamiltonian Q with orthogo- 
nal symmetry (3=1. The inset of Fig. ^ displays our re- 
sults for Ri i (s) at different b values for several system sizes: 
b = 0.02, L = 4000 (circles), b = 0.05, L = 6000 (squares), 
b = 0.08, L = 1000 (diamonds), b = 0.1, L = 4000 (up 
triangles), and b = 0.12, L = 1000 (left triangles). If the hor- 
izontal axis is rescaled by a factor 1 /b, then all data should 
collapse onto a single curve. The main panel of Fig. ^shows 
i?i i i(s) as a function of the rescaled variable s/b, thus, con- 
firming the s/b dependence of i?i i(s). The best fit of this data 
settoEq. @ gives the fitting parameter ay = 0.502±0.003, 
which is smaller than the predicted value 1 / y/ir = 0.564. Ref- 
erence 35 reported numerical values of i?x i(s) at b = 0.1 for 
two small system sizes L = 256 and 512, and, as in our cal- 
culations, their results were also in relative disagreement with 
the value l/v^ m Eq. @. The solid and dashed lines repre- 
sent Eq. @, with the fitting parameter ai,i = 0.502 and the 
predicted ai j = 1/ x/tt, respectively. Note that in Ref. l35l a 
different normalization was used in Eq. Q. 

Next we consider the d = 2 and 3 cases for which, as men- 



tioned in the Introduction, there are no analytical predictions. 
The results for Rd,i( s ) &e shown in Fig. |2]in which we were 
able to collapse both sets of data onto single curves by rescal- 
ing the normalized spacing s to the coupling constant l/b d of 
the model. Given the similarity of these results with those for 
the ID model, a curve of the form 

i^s) = 1 + S(s) - erfc (a rf)1 ^ (8) 

was fitted to the data points in this graph and found that 
a 2 ,i = 0.308 ± 0.001 and a 3 ,i = 0.208 ± 0.001. These 
fits are represented as solid lines in Fig. |2] Thus, the sys- 
tem dimensionality d of the TLCF enters via the inverse- 
coupling constant b~ d . We stress that Eq. l|8} gives a fairly 
good fit to the data. The reported data correspond to b 2 = 
0.02, L 2 = 40 2 (circles), b 2 = 0.05, L 2 = 20 2 (squares), 
b 2 = 0.08, L 2 = 60 2 (diamonds), b 2 = 0.1, L 2 = 100 2 
(up triangles), b 2 = 0.12, L 2 = 30 2 (left triangles), b 3 
0.02. /.'' = 12 3 (circles), b 3 = 0.05, L 3 = 14 3 (squares), 
b 3 = 0.08, L 3 = 12 3 (diamonds), b 3 = 0.1, L 3 = 8 3 (up 
triangles), and b 3 = 0.12, L 3 = 10 3 (left triangles). 

These results allow us to generalize the ID analytical result 
[Eq. 0] to the 2D and 3D models by simply replacing the 
inverse-coupling constant b of the Id case by the correspond- 
ing to the d-dimensional case b d . 

The observed b d dependence of Rd,i(s) in Eq. © is not 
surprising since other critical properties, such as the corre- 
lation dimension di in ID and 2D present a similar behav- 
ior toward b d . Specifically, d-2 = 1 — 1/7T& (b 3> 1) and 
d-2 = 2b (b <C 1) were derived in Refs. l33l and l34l whereas 
d 2 = 2 - a 2 /b 2 (b 2 > 1) and d 2 = c 2 b 2 (b 2 « 1) were 
numerically found in Ref. |6jJ 
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FIG. 3: As in Fig. 1, for the ID model Q with unitary symme- 
try (fS = 2). Data correspond to b = 0.02, L = 500 (circles), 
b = 0.05, L = 1000 (squares), b = 0.08, L = 2000 (diamonds), 
b = 0.1, L = 500 (up triangles), and b = 0.12, L = 1000 (left tri- 
angles). The solid and dashed lines represent Eq. J5J, with the fitting 
parameter oi,a = 0.495 and the predicted 01,2 = 2/7T, respectively. 
The inset shows the same data on the scale s. 



B. Unitary symmetry 




2 4 6 8 10 
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FIG. 4: As in Fig. 2, for the 2D (upper curve) and 3D (lower curve) 
models with unitary symmetry (J3 = 2). Data correspond to b 2 = 
0.02, L 2 = 40 2 (circles), b 2 = 0.05, L 2 = 20 2 (squares), b 2 = 
0.08, L 2 = 60 2 (diamonds), b 2 = 0.1, L 2 = 80 2 (up triangles), 
b 2 = 0.12, L 2 = 30 2 (left triangles), b 3 = 0.03, L 3 = 12 3 (circles), 
b 3 = 0.05, L 3 = 14 3 (squares), 6 3 = 0.08, L 3 = 10 3 (diamonds), 
b 3 = 0.1, L 3 = 10 3 (up triangles), and b 3 = 0.12, L 3 = 8 3 (left 
triangles). The solid lines are fits to Eq. {8}, with the fitting parameter 
fZ2,2 = 0.181 and 113,2 = 0.083, respectively. 



As in Sec. IIII Al we first analyze the ID case in order to 
compare the numerical data with the analytical result [Eq. 
@]. Figure[3]shows Ri^(s) as a function of the rescaled vari- 
able s jb. As expected, all data points collapse onto the same 
curve. The inset displays the same data on the s scale. The 
values of b and L reported are: b = 0.02, L = 500 (circles), 
b = 0.05, L = 1000 (squares), b = 0.08, L = 2000 (dia- 
monds), b = 0.1, L = 500 (up triangles), and b = 0.12, L 
1 000 (left triangles). Fitting these data to Eq. Q gives 
ax,2 = 0.495 ± 0.005, which again is small that the predicted 
value 2/7T = 0.637. To our knowledge, this is the first nu- 
merical confirmation of Eq. (|5}. The solid and dashed lines 
represents Eq. 0, with the fitting parameter 01,2 = 0.495 
and the predicted ai 2 = 2/tt, respectively. 

The results for Rd.2(s) in the d = 2 and 3 models are shown 
in Fig. |3 in which one can clearly appreciate the collapse of 
both sets of data when represented as a function of the rescaled 
variable s/b d . As in the cases of /3 = 1, a curve of the form 

RdA s ) = 1 + - ex P (-Odfi ps) (9) 

was fitted to the data points in this graph, where is 
a fitting parameter. The data reported correspond to b 2 = 
0.02, L 2 = 40 2 (circles), b 2 = 0.05, L 2 = 20 2 (squares), 
b 2 = 0.08, L 2 = 60 2 (diamonds), b 2 = 0.1, L 2 = 80 2 
(up triangles), fo 2 = 0.12, L 2 = 30 2 (left triangles), 6 3 = 
0.03, L 3 = 12 3 (circles), b 3 = 0.05, L 3 = 14 3 (squares), 
b 3 = 0.08, L 3 = 10 3 (diamonds), 6 3 = 0.1, L 3 = 10 3 (up 
triangles), and b 3 = 0.12, L 3 = 8 3 (left triangles). The 
solid lines in Fig. 0]are fits to Eq. (|9} with fitting parame- 
ters a 2! 2 = 0.181 ± 0.001 and 03,2 = 0.083 ± 0.001. Note 
that this equation gives a fairly good fit to the data. 



For the system sizes considered (L d ^> 1) we have checked 
that Rd,p{s) is an L-independent universal scale-invariant 
function, thus confirming the existence of a critical distribu- 
tion exactly at the transition. Furthermore, using the nearest- 
level distribution p(0; s), we verified that the normalized near- 
est level variances are indeed scale invariant at each critical 
point studied^ 

Before concluding the section, we summarize in TableUthe 
nonuniversal constants a,d,/3 found. Note that each value is 
different for every d and (3, thus reflecting its dependence on 
the Hamiltonian symmetry and dimensionality. The values in 
brackets were obtained from the spectral compressibility (see 

Sec. Eg. 



IV. LEVEL NUMBER VARIANCE 

This section is devoted to the calculation of the LNV 
E 2 /3 ((n)) and the spectral compressibility \ of models Q 
with the orthogonal or unitary symmetry in the strong cou- 
pling regimen, b d <C 1. Our results for \ are compared with 
the analytical prediction of Ref.l63i 

It is well known that the statistical properties of spec- 
tra of disordered one-electron systems are closely related 



TABLE I: Nonuniversal constants a^,^ of the TLCF The values in 
brackets were obtained from the spectral compressibility. 





d= 1 


2 


3 


(3 = 1 


0.502 (0.497) 


0.308 (0.291) 


0.208 (0.232) 


2 


0.495 (0.483) 


0.181 (0.206) 


0.083 (0.105) 



5 



to the localization properties of the corresponding wave 
functions>A2ii2i£ In particular, Ref. |63l based on a Brown- 
ian motion level approximation combined with an assumption 
concerning the decoupling of the energy levels and eigenfunc- 
tion correlations, derived the following relation for multifrac- 
tal eigenstates 



X 



2 I d 



(10) 



where d 2 is the correlation dimension of the critical wave 
functions. According to this result, the spectral compressibil- 
ity x should tend to 1/2 at the limit of very sparse eigenstates 
d 2 — > 0, and not to the Poisson value \ = 1, Although 
Eq. J10I has been widely confirmed at weak multifractality 
(weak coupling limit) its validity at strong coupling (small 
d 2 ) raises many doubts. More precisely, in Ref. |35j it was 
analytically shown that at the limit of small 6, % — > 1 for 
the one-dimensional model Q. This tendency has been also 
numerically demonstrated in Refs. 35 and 64, and for the An- 
derson model in d > 

We are therefore especially interested in the calculation of 
X for d — 2 and 3 in order to check whether or not Eq. dlOi 
adequately describes its behavior at strong coupling b d <C 1. 
The compressibility x can be expressed through the TLCF as 



X 



ds [Rd,p{s) - 1] 



(11) 



Substitution of Eqs. (@), (0, ©, and l|9} into Eq. (JTT) yields 
the spectral compressibility for orthogonal ((3 = 1) as well as 
for unitary symmetry (j3 = 2) 



X = 1 - Cd,f> b d 



(12) 



where Cd,i — Ijad^wpK and Cd. 2 — \J ^/o-d.2- Equation (I12> 
constitutes a generalization to d dimensions of the ID ana- 
lytical estimates of Ref. |35j Notice that at the limit of very 
strong coupling b d — > 0, \ tends to the Poisson value 1 in 
disagreement with Eq. JlOi . 

An alternative way to directly calculate \ is from the 
asymptotic behavior of the LNV [Eq. 0]. We will show 
that for b d <C 1, relation H2\ is satisfied, thus, giving further 
support to the proposed relations for the TLCF [Eqs. (|8j and 

In Fig. |5] we show the computed £ 2 1 (("■)) f° r me 2D or- 
thogonal ensemble (J3 = 1) at different disorders and system 
sizes: b 2 = 0.02, 1? = 40 2 (right triangles), b 2 = 0.02, L 2 = 
60 2 (stars), b 2 = 0.04, L 2 = 40 2 (down triangles), b 2 = 
0.04, L 2 = 60 2 (diamonds), b 2 = 0.08, L 2 = 20 2 (left tri- 
angles), b 2 = 0.08, L 2 = 60 2 (squares), b 2 = 0.1, L 2 = 20 2 
(up triangles), and b 2 = 0.1, L 2 = 100 2 (circles). Note that, 
for each value of b 2 , S| 1 ((n)) is L independent, which is a 
sign of criticality. There is a clear gradual tendency in the 
large (n) region of £?, x((n)) toward the Poisson limiting re- 
sult (dashed line) as the inverse coupling constant b 2 of the 
model is decreased. We have checked that for (3 = 2 and for 
the 3D case the behavior is quite similar. The straight lines, 
whose slopes correspond to the values of \, are fits to Eq. Q 
with the fitting parameters summarized in Fig. [6] 




FIG. 5: Level number variance E|,i {{n}) for the 2D orthogonal 
ensemble (/3 = 1) at different disorders and system sizes: b 2 = 
0.02, L 2 = 40 2 (right triangles), b 2 = 0.02, L 2 = 60 2 (stars), 
b 2 = 0.04, L 2 = 40 2 (down triangles), b 2 = 0.04, L 2 = 60 2 (dia- 
monds), b 2 = 0.08, L 2 = 20 2 (left triangles), b 2 = 0.08, L 2 = 60 2 
(squares), b 2 = 0.1, L 2 = 20 2 (up triangles), and b 2 = 0.1, L 2 = 
100 2 (circles). Solid lines are fits to Eq. jTJ, and the dashed line 
corresponds to the Poisson result. 



The b d dependence of the spectral compressibility x, as ob- 
tained from the previous fits, for the ID (squares), 2D (cir- 
cles), and 3D (diamonds) model Q with orthogonal (solid 
symbols) or unitary symmetry (open symbols) is depicted in 
Fig. [6] This clearly shows that for almost all values of b d 
reported, x I s greater than the maximum value of 0.5 pre- 
dicted by Eq. dlOl i. The solid lines are fits to the form 
X = Xo~ c d./3b d with the fitting parameters ci.i = 2.27±0.14, 
ci.2 = 2.55 ± 0.13, c 2 .i = 3.88 ± 0.12, c 2 , 2 = 3.90 ± 0.26, 
c 3 'i = 4.87 ± 0.19, and c 3 , 2 = 5.46 ± 0.20. Using these 
values of Cd,p and taking into account their relation with the 
nonuniversal constants dd,p of the TLCF, Cd.i = 2/ad.iV^ 
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FIG. 6: b dependence of the spectral compressibility x for the ID 
(squares), 2D (circles), and 3D (diamonds) models Q with orthogo- 
nal (solid symbols) or unitary symmetry (open symbols). Solid lines 
are fits to the form \ = Xo — Cd,/3b d . 
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and Cd,2 — \/ft/ad.2, we can obtain ad,p from a different 
quantity. The corresponding values are given in brackets in 
TablelU in which one can appreciate the good agreement with 
those obtained from the TLCF. Based on these results, one 
can say that, in the case of strong coupling, the Brownian mo- 
tion level approximation breaks down and Eq. dlOt becomes 
invalid for all values of d considered. 



V. SUMMARY 

We present the first numerical results for the two-level 
correlation function Rd,p(s) of noninteracting electrons on 
a d-dimensional disordered system with long-range transfer 
terms. Models with orthogonal or unitary symmetry at small 
values of the inverse-coupling constant b d have been consid- 
ered. The ID analytical results [Eqs. @ and Q] are con- 
firmed (except for the numerical constants). We also found 



that the ID formulas are valid for the 2D and 3D models if the 
inverse-coupling constant b is replaced by the corresponding 
to the d-dimensional case b d . Another important result pro- 
vided by our numerical calculations is that the spectral com- 
pressibility, which is found to be close to 1 in the limit b d — ► 0, 
does not satisfy the relation JlOi . The proposed Eqs. (|8j 
and (|9) are based on numerical results and, at present, should 
be considered as conjectural. So, further analytical work is 
needed to check these forms of the TLCF and their origin in 
the model Q. 
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